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For  a  dispersive,  non-magnetic ,  electrically  conducting  medium,  Maxwell's 
equations  for  the  propagation  of  electromagnetic  waves  are  given  in  terms  of 
dimensionless  variables  by 


VxE  -  -  H.  , 

VxH  -  D  , 

(1) 

-  ~t 

D  =*  (E  +  P) 

—  — t 

(2) 

The  vectors  D(x,t),  E(x,t),  H(x,t)  and  P(x,t)  are  proportional  to  the  electric 

displacement,  the  electric  field  strength,  the  magnetic  field  strength,  and 
the  polarization  vector,  respectively.  The  constitutive  properties  of  the 
medium,  which  are  experimentally  determined  give  the  polarization  vector  in 
(2)  as  a  linear  functional  of  E.  A  general  expression  is  given  by  the 


differential  constitutive  law. 


P(x,t)  -  2  P  (x  t) 

j-1  ~J 


S  or  .  “  —  S  0.  ~ Z  ,  J-1.2....N  .  (3b) 

n=0  nJ  3tn  m-0  mJ  3tra 

The  specified  constitutive  coefficients,  and  ^ ,  m  -  0,1,... s, 

n  =  0,1 . r  and  j  =  1,2,...,N  are  experimentally  determined. 

More  generally,  we  express  P  as  a  linear  functional  of  E  by  the  heredity 

integral , 


V  *  ^ 


P(x,t)  =^E(x,t)  *  J  E(x,t')€(t-t')dt'  , 


where  the  "memory"  function,  e(r),  which  is  experimentally  determined, 


vanishes  for  r  <  0. 


We  have  first  analyzed  the  propagation  of  finite  jump  discontinuous 


solutions  of  the  system  (1) - (3)  using  both  the  method  of  weak  solutions  [1,2] on/ _ 

tty  Codes 

and  progressing  wave  expansions  [1,3]  for  the  Debye  [4]  model  of  the  ,* and/or 

uist  spooial 

U  I 


dispersive  medium.  This  model  is  believed  to  be  a  good  approximation  of  the 
dispersive  properties  of  biological  materials.  Thus,  each  satisfies  a 

first  order  ordinary  differential  equation 

fc  E]  +  “j?j  -tjS  ■  J  -  i.2 . »  •  <4> 

The  "constitutive"  coefficients  a.  ,a_ . a  and  fl,  ,j9  are  to  be 

12  n  1^2  n 

determined  experimentally,  or  from  quantum  mechanical  calculations.  Then, 


pt  -a . 

P.  -  J  E(x,t')^e  J 


a.(t-t') 


t  N  -a.(t-t') 

and  P  -  j"E(x,t)  £  fi.e  J  dt' ,  so  that  the  memory  function  for  the  Debye 

~  j-1  J 

model  is  given  by  the  linear  combination  of  exponentially  decaying  functions, 

N  -a .  r 
e(r)  -  £  fi.e  J  . 

j-1  J 

The  progressing  wave  expansion  for  the  solution  of  (1) , (2)  and  (5)  is  a 
representation  in  the  form 

CO  00 

E-  £  f  [*(x,t)]Em(x,t)  ,  H  -  £  f  [*(x,t)]Hm(x,t)  , 

m-0  m  ~  “  ~  ra-0  m  ~  ~ 


P.  -  £  f  [*(x,t)]P?<xft)  , 

~J  m-0  m  ~  “ 

where  £q(.4>)  is  an  arbitrary  function  of  the  single  variable  ^(x,t)  and 
f^,f2 .  are  required  to  satisfy, 

£'  -  f  ,  ,  m  -  1,2 .  (7) 

m  m-1 

Then  we  insert  (6)  into  (1),  (2)  and  (4)  and  use  (7)  for  the  derivatives  of 
f  .  By  equating  to  zero  the  coefficient  of  each  function  fm(^) .  we  obtain  a 


sequence  of  equations  to  determine  <f>,  Era,  Hra  and  P™.  For  example,  if  we  set 
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<^(x.t)  —  t  -  V*(x)  »  then  satisfies  the  eikonal  equation  of  geometrical 
optics , 


iw-r  - 1  , 


(8) 


and  the  coefficients  E^,  E^",  etc.,  satisfy  transport  equations,  which  reduce 


to  first  order  partial  differential  equations  along  the  characteristics  of 
(8).  Solving  the  transport  equations  explicitly,  shows  that  E^  and  decay 


exponentially  along  the  rays. 

The  choice  of  the  waveform  fg(<£)  Is  determined  by  the  type  of 

discontinuity  to  be  analyzed,  and  the  form  of  the  pulse.  For  example,  to 
recover  the  results  of  our  previous  analysis  of  the  propagation  of  jump 
discontinuities,  we  select  f^  to  be  the  Heaviside  step  function  i.e. 


f0  -  H(t-tf) 


For  the  Lorentz  model  [4]  of  dispersion  the  polarization  vector  is  given 
by  (3a)  where  the  P.  now  satisfy 


32P. 


3P. 


+  a,  TT-  +  7 i  P* 


J  ~3 


P  E 


(ID 


and  the  constitutive  coefficients  a  ,  and  are  to  be  experimentally 

determined.  By  solving  the  transport  equations  for  the  amplitude  of  the 
progressing  wave  expansion  of  Maxwell's  equations  using  the  Lorentz  model,  we 
find  that  the  amplitudes  do  not  decay  exponentially  along  the  ray.  Thus,  we 
are  now  studying  more  general  models  of  the  dispersive  media,  e.g.  (3),  to 
determine  the  qualitative  features  for  which  the  amplitudes  do,  or  do  not 
decay  exponentially.  Preliminary  results  suggest  that  a  classification  of 
dispersive  media  is  obtained  depending  on  the  relative  orders  of  the 
differential  operators  on  the  right  and  left  sides  of  (3.3b). 
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Finally,  using  the  progressing  wave  expansion,  we  have  determined  the 
scattering  of  pulses  from  dispersive  half  spaces  and  from  dispersive  targets. 

A  paper  entitled  "Progressing  Wave  Expansions  for  Temporally  Dispersive 
Electromagnetic  Waves  and  the  Classification  of  Dispersive  Media"  is  currently 
in  preparation. 
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